Non-classicality thresholds for multiqubit states - numerical analysis 
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States that strongly violate Bell's inequalities are required in many quantum-informational pro- 
tocols as, for example, in cryptography, secret sharing, and the reduction of communication com- 
plexity. We investigate families of such states with a numerical method which allows us to reveal 
non-classicality even without direct knowledge of Bell's inequalities for the given problem. An 
extensive set of numerical results is presented and discussed. 



I. INTRODUCTION 

Until the discovery of the protocol of entanglement- 
based quantum cryptography [1[ , Bell's theorem was rel- 
evant only in discussions concerning foundations of quan- 
tum mechanics. It established means to evaluate the 
most reasonable class of hidden variable theories, i.e., lo- 
cal and realistic theories. Once the theorem proved useful 
in quantum cryptography and other quantum communi- 
cation schemes (secret sharing [2| , reduction of commu- 
nication complexity [3-13: etc.), it gained an additional 
technical status. Today, we are not only interested in 
showing that given quantum correlations do not admit a 
local hidden variable model, but also, with Bell inequal- 
ities, we can now prove the usefulness of states in quan- 
tum protocols. The strength of such violations indicates 
resistance to noise or to decoherence. 

However, only in simple cases, like two dichotomic 
measurements per observer, one can pinpoint the full 
set of such Bell-type inequalities, and put them into a 
manageable form of just one single non-linear inequality 
0| • I'^ other cases we know usually only some subsets 
of tight inequalities for the given problem. These subsets 
still contain a vast number of different inequalities. 

One of the methods to avoid the problem of both miss- 
ing inequalities and handling the exploding numbers of 
known ones is to address the problem numerically, with- 
out a direct use of Bell's inequalities |8j. It is very easy 
to notice that the problem of existence of a local hidden 
variable model for the given set of quantum probabilities 
is a form of a linear programming problem (see Sec. Ill| . 
Thanks to such an approach it was possible to show for 
example that non-classicality of two-qudit correlations is 
higher than for qubits, and that it increases with the 
dimension Q- 
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Below we present how linear programming can be har- 
nessed to study non-classicality of quantum correlations. 
This approach was put in the form of a computer code: 
Steam-roller. We will present the numerical results 
in Sees. Ill and IV. Here we particularly aimed at 
multi-setting scenarios for multi-qubit states which are 
known to possess interesting kinds of non-classical be- 
havior. Most of our results cover regions for which the 
full set of Bell's inequalities is unknown. Even the results 
for two settings per observer cover a new territory: they 
are produced for the full set of probabilities for the given 
iV-qubit problem. In the case of the WWWZB inequal- 
ities [y, 0, [i3| only highest order correlation functions 
[111 are taken into account. These do not reflect inter- 
ference effects that may occur between subsets of par- 
ticles. Thus the WWWZB set, despite being complete 
for highest order correlation functions (in the case of a 
two setting choice per observer), is not complete for the 
full description of quantum phenomena. This full set of 
probabilities for the given problem gives full knowledge 
about an experiment, and therefore contains all available 
information. 



II. LINEAR PROGRAMMING 

In the most general Bell's experiment, N observers per- 
form measurements on a given state pN- Each observer 
can choose between rrii arbitrary dichotomic observables 
(i = 1,...,A^). Their measurement results will exhibit 
correlations, not describable by classical statistics (i.e., 
by classical realistic models). For brevity let us consider 
a case with two observers (Alice and Bob) and two mea- 
surement settings per side. In this case the local realistic 
models are equivalent to the existence of a joint probabil- 
ity distribution pir{ai,a2,hi,h2), where a^ = ±1 denotes 
the result of the measurement of Alice's i-th observable 
(Bob's results are denoted by bk — ±1). Quantum pre- 
dictions for the probabilities are the marginal sums of the 



above joint probability distribution pir 

PQM{ra,rb\Ai,B2)= ^ Plr{r a, 0.2, hi, n) 

a2,bl— ±1 

PQM{ra,n\A2,Bi) ^ ^ Piriai,ra,n,b2) 

ai,f)2=±l 

PQMira,n\A2,B2) ^ ^ p;r(ai,ra,6i,r6), (1) 
ai,bi— ±1 

where PgAfC'^a, '"fcl^i, S^) denotes the probability of ob- 
taining the result r^ by Alice and ri, by Bob, when they 
measure the observables Ai and Bk, respectively. It 
has been shown that for some entangled states, no local 
realistic probability distribution p;^ (01,02, fei, 62) exists, 
which could satisfy the set of equalities (H]). Therefore, 
no local realistic model can reproduce the predictions 
of quantum mechanics for these entangled states. This 
statement is known as Bell's theorem [13 |. 

This is one of the ways one can express Bell's theorem. 
The more frequent approach is to find certain (Bell) in- 
equalities which are satisfied by local realistic models, 
but violated by some quantum predictions. More im- 
portantly, there are many statements on the "strength" 
of a violation of such inequalities which use such terms 
as "amount of violation" or "factor of violation." How- 
ever, such measures cannot be used for direct compar- 
ison between different Bell inequalities or with a differ- 
ent method, such as the one used in our paper. Thus, 
following Q, we shall use a "violation" parameter Vcrit 
which seems to be much more objective: we shall seek 
the amount of random ( "white" ) noise admixture that is 
required to completely hide the non-classical character of 
the original correlations for the given state. 

If we mix some amount of white noise to the two-qubit 
state p, we obtain a state described by the following den- 
sity operator: 



p{v) = vp - 



1-v 



22 



-1 



82 



(2) 



The quantmn probability PQj^j{ra,rb\Ai, Bk) — 
P{ra,rh\Ai,Bk) for the state p{v) reads 

P{ra,n\A,,Bk)^vPQM{ra,n\A„Bk) + ^il-v). (3) 

The parameter v is the visibility of the state, and obvi- 
ously (1 — v) is the amount of noise admixture. For v = 
the equalities ([TJ are obviously satisfied (white noise has 
a local realistic model). For v = 1 the equalities ([T]) can 
not be satisfied for some entangled states p (and for a par- 
ticular choice of observables) . For such states there exists 
the critical visibility Vcrit that for v ^ Vcrit there exists 
a local realistic probability distribution pir{ai, 02, 61, 62) 
that satisfies the set of equalities ([T]). 



A set of 16 probabilities PQM{ra,rb\Ai,Bk) corre- 
sponds to a single point in a 16-dimensional space. Those 
sets which can be described using local realistic mod- 
els form a convex geometric figure called a Bell-Pitovsky 
polytope in this 16-dimensional space. The facets of the 
polytope are equivalent to tight Bell inequalities. For 
V — Vcrit the local realistic probabilities pirifli, 02, &i, 62) 
lead to a set of probabilities PQM{Ta,T'b\Ai,Bk) which 
form the coordinates of a geometric point which lies in 
a facet of the polytope (i.e., they saturate a Bell in- 
equality). It has to be noted, however, that our nu- 
merical method does not lead to uncovering particular 
Bell inequalities which are violated by an investigated 
state. We think that, the power of the method lies in the 
fact that one docs not need any knowledge at all of the 
forms of Bell inequalities when analyzing non-classicality 
of quantum-entangled states. Seeking violated Bell in- 
equalities is highly inefficient, since even in the simplest 
cases one ends up with an exploding number of Bell in- 
equalities, the bulk of them trivial. We skip this cumber- 
some step entirely, see the following paragraphs. 

As mentioned, critical visibility Vcrit depends on the 
particular set of observables Alice and Bob choose from. 
We can parametrize any dichotomic observable Xi by the 
two angles Of and (j)f {X ~ A, B) in the following way: 



X = 



)x{ 



)x{ 



\±)x =cos(±7r/4-f 6if)|0)x + e''^- sin(±7r/4 + 6if )|l)x. 

(4) 

Now, the quantum probability can be calculated: 
PQM{ra,rb\Ai,Bk) = Tr{p\ra)A{ra\ <E) \rb)B{rb\)- Hence 
we can define the critical visibility function of the angles 

Vc„t{Ot,'l^t,Otcl^tOf,cl,f,ei,<j>^) e [0,1]. (5) 

Our task is to find, for a given state p, the minimal 
critical visibility. If the two observers A and B perform 
two measurements each, the probabilities ([3]) are then 
parametrized by the angles: 

P{ra,rb\A,Bk) = P{ra,rb,9f,cf>f,9^,^^). (6) 

The critical visibility is obtained by maximization of the 
visibility until the set of equalities ([Ij can no longer be 
satisfied. This is done by means of linear programming. 
We adopt 

Po =Pir{-, ■■■,-,-), 
Pi =pir{-, ■■■,-,+), 



Pn-l =Plri+, ■■■,+,+), 



(7) 



where 



omi-t-. 



.+mjv g^j-^jj hereby give a complete de- 



scription of the linear programming problem to be solved 
(based on Eqns. ([J) and ([3])) for the case with Alice and 
Bob: 



y^ Plr{ra,0,2.n,b2 



02,62 — =tl 



l-v 



+ i;P(ra,r,,0f,0f,0f,0f) 



a2,&i=±l 



1-V 



Y^ Pir{ra, a2,bi,n) = ^^ + vP{ra, n, et (f>t Oi,(t>i), 



y^ Pir{ai,ra,rb,b2) 



0.1 ,62— ±1 



l-v 



l-v 



+ vPira,n,9^,(bt,0r,(bf) 



Y Pir{ai,ra,h,n) = —2- +wP(ra,rf„6'^,(/)^,6'f ,(/)f), 



ai ,61— ±1 
n-1 



^Pfc = l; O^ws;!; < pfc «: 1, 

/c=0 

z{pQ,Pi,...,Pn-l,v) ^V, 



(8) 



where z is the (trivial) function to be maximized. 

In the general case the four expressions above evalu- 
ate to c = 2^mi • • • 771 jv equalities (e.g. for 3 qubits and 
6 measurement settings per side, there are 1728 equal- 
ities). The above is known as the canonical form of 
the linear programming problem. We use the revised 
simplex method [13l | from the GNU Linear Programming 
Kit [ij] to solve this problem. The top constraints form 
a simplex in the (n + l)-dimensional space spanned by 
Po, • • • ,Pn-i and V. The algorithm traverses the vertices 
of the simplex until it finds the optimal solution. The 
returned solution is the global maximum of z - the crit- 
ical visibility we seek. The minimization over the angles 
ef, 4't,0f,(l)f,e^^, e^,<j)^ is realized by the downhiU 
simplex method J15l | from the SciPy package [3 • 

Please note that the complexity of the linear program- 
ming problem is defined by n and c ~ values which in- 
crease exponentially with the number of observers N or 
the number of settings rrii per observer. We are, there- 
fore, limited to the computational capabilities of today's 
numerical machines when solving this problem. We used 
a machine with an Intel Pentium D CPU 3.20GHz pro- 
cessor for our computations. Some of the results in this 
paper took several weeks to compute. Computations of 
results for bigger problems can be anticipated to com- 
plete in a reasonable time when using faster machines in 
the future. 

One might wonder about the credibility of the results 
presented in this paper. There are two aspects which 
raise doubts, namely, numerical precision and local min- 
imum risk when computing the minimum of the criti- 
cal visibility function. Please note that there is no local 
maximum risk when computing the value of the critical 
visibility function - the simplex algorithm guarantees it 
will always find the global maximum as described in [13j . 
As for the precision of the results, we will not conduct an 
extensive analysis but will plainly state that this is not 
an issue. For all the results which can be calculated ana- 
lytically, Steam-roller generated results are compliant 
on all 14 displayed significant digits. For example, for a 



I 

two-qubit Greenberger-Horne-Zeilinger (GHZ) state, the 

analytically-derived minimal critical visibility is -^ and 
Steam-roller gives 0.70710678118655. One of the rea- 
sons for such outstanding precision is that the numeri- 
cal errors do not propagate between consecutive critical 
visibility function invocations. Another is that the linear 
programming library used properly handles numerical in- 
stabilities. The local minimum risk when computing the 
critical visibility function is more of an issue. There is 
no guarantee we have reached the global minimum, and 
Steam-roller did fall into a local minimum from time 
to time. One should, therefore, treat all the results in 
this paper in the following way; it is guaranteed that a 
given state has a minimal critical visibility not greater 
than the one presented. As for the likelihood that it is 
also not smaller, we can only state that there's a good 
chance it is because Steam-roller provided us with a 
correct result in 99 cases out of 100. 



III. BOUNDS ON LOCAL REALISM 

We applied the numerical method to different types of 
quantum states, which are often discussed in the context 
of quantum information studies, namely, 

• the generahzed GHZ state [13, [II: 
|GHZ(a))7v = cosa|0 • • • Q)n + sinall • • • 1)n, 

• the W state [H: 

|W)w - ^(|10 • • • 0)w + |010 • • •0)jv + - • - + 10 • • -01)^), 



I the four-qubit singlet state [3, [20| : 
1*4) = ^(10011) + 11100)) 

(10101) + 10110) + |1001) + |1010)), 



• the six-qubit singlet state [20|, |21| : 

|*6> = ^\GHZ-)e + 1{\W)3\W)3 - \W)3\W)3), 

where \GHZ-)e = ^(|000111) - |111000)), and 
\W)3 is the spin- flipped \W)3, 

• the symmetric Dicke state |22| : 



\D 



{N/2) 

N 



N 

N/2 



-1/2 



Y^ I0---01---10---0) 

N/2 



N, 



permutations 



the four-qubit cluster state [23|: 



1 



ICluster) = -(10000) + |0011) + |1100) - |1111)). 



We calculated the critical visibility for an increasing 
number of different settings per side. All results are pre- 
sented in Tab. HI They lead to the following observations. 

Observation 1. An explicit form of tight Bell inequal- 
ities for probabilities is known only for the following 
experimental situations: 2x2 [13], 3x3 [H, HI], 2x2x2 
[25l |26|, whereby ixj means i(j) alternative measure- 
ments on Alice's (Bob's) side. In these cases the nu- 
merical method gives the same critical visibility as the 
analytical expressions. For states (#1 - #5) (see Tab. H]) 
we observe such a correspondence. 

Observation 2. Local realistic correlations in the N- 
qubit experiments with two alternative measurement set- 
tings per side are bounded by WWWZB inequalities 
0)0) ll3- I^ t^is case, the numerical method should give 
better or at least the same results. We obtain the same 
results for the following states: GHZ(7r/4)Ar (#1, #2, 

#13, #31, #35, #39), *4 (#21), i?f ^ (#25), Cluster4 
(#29), ^6 (#37), Of^ (#38). Stronger constraints on 
local realism than predicted by the WWWZB inequalities 
are observed for the following states: 

• |W)Ar (#17, #33, #36, #40). The corresponding 
critical visibilities obtained by means of WWWZB 
inequalities are equal to: 



0.6325, 



- 0.6202, i;-f 



'Wz 



= 0.6565, z;^f 



= 0.6141, and 



0.6089, and are higher than obtained by the numer- 
ical method by about 2% (for A^ = 3) to 208% (for 
N = 7). The values (#5, #17, #33) were previ- 
ously published in f27'|. 

• GHZ(7r/12)Ar (#3), GHZ(7r/180)Ar (#4). These 
states do not violate WWWZB inequalities at all 

[11 [H. 

One of the simplest examples showing the advantage of 
Bell inequalities based on probabilities over correlation 
functions is a three qubit state, which is the product of 
the two qubit singlet state and white noise \'ip~)ab{4'~\® 
Ic- The state does not violate the WWWZB inequalities 



No. 


iV State 


mi X 7712 X • ■ • 


X 777 JV 


Vcrit 


#1 


2 |GHZ(7r/4))2 


2x2 - 10x10 




0.7071 


#2 


3 |GHZ(7r/4))3 


2x2x2 - 5x5x5 




0.5000 


#3 


|GHZ(7r/12))3 


2x2x2 - 4x4x4 




0.8165 


#4 


|GHZ(7r/180))3 


2x2x2 - 3x3x3 




0.9988 



#5 

#6 

#7 

#8 

#9 

#10 

#11 

#12 



#25 
#26 
#27 
#28 



\W)3 



2x2x2 
3x3x2 

4x4x2 
5x5x2 
6x6x2 
3x3x3 
4x4x4 
5x5x5 



\D':') 



2x2x2x2 
3x3x3x2 
4x4x4x2 

3x3x3x3 



0.6442 

0.6330 

0.6240 

0.6236 

0.6236 

0.6048 

0.6013 

0.6007 



#13 4 


|GHZ(7r/4))4 


2x2x2x2 


0.3536 


#14 




3x3x3x2 


0.3536 


#15 




4x4x4x2 


0.3536 


#16 




3x3x3x3 


0.3408 


#17 


\W)i 


2x2x2x2 


0.5469 


#18 




3x3x3x2 


0.5011 


#19 




4x4x4x2 


0.4923 


#20 




3x3x3x3 


0.4900 


#21 


l*4> 


2x2x2x2 


0.5303 


#22 




3x3x3x2 


0.5012 


#23 




4x4x4x2 


0.4948 


#24 




3x3x3x3 


0.4823 



0.4714 
0.4646 
0.4630 

0.4393 



#29 
#30 


Cluster) 4 


2x2x2x2 
3x3x3x3 


0.5000 
0.4472 


#31 5 
#32 


|GHZ(7r/4))5 


CO to 
X X 
CO to 
X X 
CO to 
X X 
CO to 
X X 
CO to 


0.2500 
0.2280 


#33 
#34 


\W), 


CO to 

X X 
CO to 
X X 
CO to 
X X 
CO to 
X X 
CO to 


0.4300 
0.3462 


#33 
#34 


Cluster) 5 


CO to 
X X 
CO to 
X X 
CO to 
X X 
CO to 
X X 
CO to 


0.3333 
0.3333 


#35 6 


|GHZ(7r/4))6 


2x2x2x2x2x2 


0.1768 


#36 


\W), 


2x2x2x2x2x2 


0.2927 


#37 


1*6) 


2x2x2x2x2x2 


0.3536 


#38 


l/^D 


2x2x2x2x2x2 


0.2827 


#39 7 


|GHZ(7r/4))7 


2x2x2x2x2x2x2 


0.1250 


#40 


\Wh 


2x2x2x2x2x2x2 


0.1975 



TABLE I. The critical visibilities for various types of quan- 
tum states and experimental situations (number of settings). 
If w > Vcrit, there does not exist any local realistic model 
describing quantum probabilities of experimental events. 



for three qubits, whereas the numerical method detects 
violation of local realism and gives the critical visibility 
equal to 0.7071. This number is equivalent to the visibil- 
ity, which is necessary to violate the Clauser-Horne in- 
equality by the \ip~)ab state alone. The advantage comes 
from the fact that probability methods use the whole 
knowledge obtained from an experiment. 

Observation 3. Explicit forms of tight correlation func- 
tion Bell inequalities |20| for three qubits and three set- 
tings per side were given in [30|. The numerical method 
gives the same result in (#2). In the case of (#3, #4), 
the inequalities of [30| are not violated. The critical vis- 
ibility for violation of the inequalities by the W3 state 
(#5) is equal to 0.6547 (about 2% above the visibility 
predicted by the numerical method). 

Observation 4- Explicit forms of tight correlation func- 
tion Bell inequalities for many qubits {N > 3) and many 
measurement settings (rui > 2) are still unknown. Tight 
Bell inequalities were constructed for some particular 
cases (e.g. 2^-1x2^-1x2^-2^. . .x2 [Sll) only. There 
is also the family of the correlation Bell inequalities for 
many qubits and an arbitrary number of settings [32| . 
Unfortunately, these inequalities are not tight. In this 
case, the numerical method gives the strongest known 
constraints on a local realistic description. 

• First, let us consider the 4x4x2-type inequality of 
[3l|. It is violated by the GHZ(7r/4)3 state with a 
critical visibility equal to 0.3536 which agrees with 
(#2). However, for the W3 state the inequality is 
violated if the visibility is higher than 0.6547 }31j , 
whereas the numerical method reveals the critical 
visibility to be lower by about 4.9% (#7). 

• One can also compare the numerical method with 
so called geometric inequalities [32]. For instance, 
the critical visibility of the four qubit GHZ state 
to violate the three-setting geometrical inequality 
'I'fe^m'' = 0.3421 is close to (but higher than) that 
expected by the numerical method (#16). How- 
ever, for the four-qubit W state the three-setting 
geometrical inequalities are violated with a critical 
visibility w^o„ — 0.6843, whereas the numerical 
method provides a drastically (about 28%) better 
result (#20). 

Observation 5. One can conjecture that the critical vis- 
ibility for the three-qubit W state for an arbitrary large 
number of settings is equal to 6/10 (#12). 

Observation 6. In many cases, the critical visibility for 
impossibility of a local realistic description decreases with 
the number of settings. The highest visibility drop (19%) 
is observed for the five-qubit W state (#33, #34), 10% 
for the four-qubit W state (#17, #20), and less for W3 
(#5, #10 - #12), GHZ(V4)4 (#13, #16), GHZ(7r/4)5 
(#31, #32), Cluster4 (#29, #30), ^4 (#21, #24), and 

Dl (#25, #28). Interestingly, for two- and three-qubit 
GHZ states we do not observe this effect. Acording to [33| 
it could be that the effect appears with higher numbers 
of measuremaent settings. 



Observation 7. The minimal number of measurement 
settings determines a lower bound on the critical visibil- 
ity. In (#5 - #10, #13 - #16, #17 - #20, #21 - #24, 
#25 - #28) the last observer has only two alternative 
measurement settings. Increasing the number of mea- 
surement settings at the other sides does not result in 
the decrease of the critical visibility below the value for 
three settings per each observer. 

Observation 8. The 7V-qubit GHZ(7r/4) state reveals 
the strongest non-classical properties (the lowest critical 
visibility for an arbitrary number of A'' < 8). However, 
one can observe that the difference between the critical 
visibility for the GHZ state and the W state (or the clus- 
ter state) decreases with the number of qubits. Accord- 
ing to [23, we know that the W state leads to stronger 
nonclassicality than GHZ states for N > 10. Finding an 
example for A^ < 10 is still an open problem. Due to 
computer power limitations, we were not able to recover 
the result of Sen et. al or find another one for A^ ~ 10. 



A. Bound entanglement 

We also analyze some classes of bound entangled 
states, namely. 



• the three qubit Bennett state [3J]: 
1 



PBennett 



r(i 



4 

E 



10.) (0.1), 



where |(/.i) - |O1+),|02) - |1 + O),|03) 
01),|</)4) = |-— ),|±) = (|0)±|1))/V2, 



• the A-qubit Diir state [3 
1 



Dur 

Pn = 



with 



/ 1 ^ 

\ k=l 



Pk 



-i=[|0)i...|0)Ar + e^""|l)i...|l)^], and 



V2 



Pk being a projector on the state |0)i...|l)fc...|0)Ar 
with "1" on the fcth position (i\. is obtained from 
Pk after replacing "0"s by "l"s and vice versa). 



• the generalized A^-qubit Smolin state j3 



Smolin 

Pn 



1 

2^ 






where ai are the Pauli matrices. 

The results are presented in Tab. |IT] and lead to the 
following observations. 

Observation 9. Violation of local realism by the bound 
entangled Smolin state (for four and six qubits) . The crit- 
ical visibilities obtained by the numerical method (#43, 
#47) recover the results of [33 . If we go to three alter- 
native settings per side, the critical visibility for the four 
qubit Smolin state decreases (#44). 



No. iV state 



mi X m2 X 



X mjv Vcrit 



#41 3 |Bennett)3 



2x2x2 - 5x5x5 



1.0000 



#42 4 |Dur)4 
#43 |Smolm)4 
#44 



2x2x2x2 - 3x3x3x3 1.0000 
2x2x2x2 0.7071 

3x3x3x3 0.6986 



#45 5 |Dur)5 



2x2x2x2x2 -3x3x3x3x3 1.0000 



#46 6 |Dur)6 
#47 |Smolm)6 



2x2x2x2x2x2 
2x2x2x2x2x2 



1.0000 
0.7071 



#48 7 |Diir)7 



2x2x2x2x2x2x2 



1.0000 



TABLE II. The critical visibilities for some bound entangled 
states and experimental situations (number of settings). If 
V > Vcrit, there does not exist any local reahstic model de- 
scribing quantum probabilities of experimental events. 



Observation 10. In [S^, |3^, |33| the problem of the vio- 
lation of Bell-type inequalities by the Diir state is consid- 
ered. A simple experimental situation involves 7 qubits 
and 3 settings [37i| , 6 qubits and 5 settings [32i| , or 8 qubits 
and 2 settings per side [35j]. We only verified the case of 
7 qubits (#48). Using the numerical method we can not 
find a violation of a local realistic model for two settings 
per side. It suggests that the requirement of three set- 
tings is reasonable. We also partially verified the cases 
of 4 (#42), 5 (#45) and 6 (#46) qubits. We can not find 
a local realistic model for two measurement settings per 
side (and three setting in the case of 4 and 5 qubits). The 
problem of a violation of local realism by predictions of 
the four and five qubit Diir state, and the six qubit Diir 
state with the number of measurement settings less than 
5 is still an open question. 

Exemplarily, in App. [B]we give a local realistic model 
for quantum probabilities in a Bell experiment for the 
four qubit Diir state. The observers choose between two 
measurement settings which are optimal for the case of 
the GHZ state. 



IV. FROM DETECTOR CLICKS TO DENSITY 
MATRIX 

In this section we use the numerical method for ex- 
perimentally observed states. To this end, we need their 
density matrices. Methods that are used to obtain the 
density matrix from the measured data are called quan- 
tum state tomography. Different ways for quantum state 
tomography have been developed (one of the first works 
was [38]; a variety of methods is summarized in [39|). 
From a decomposition of the density matrix into projec- 
tors, one can easily express the density matrix in terms of 
probabilities for detecting a certain coincidence. Relative 
frequencies obtained in a measurement are, however, sub- 
ject to Poissonian counting statistics. Because of these 
uncertainties, the deduced density matrices might be un- 
physical, e.g. one might find negative eigenvalues. Such 
raw data are not useful for the numerical method and 
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FIG. 1. (Color online) The bars represent for several families 
of quantum states the critical visibilities, which are necessary 
to falsify a local realistic description in two setting experi- 
ments. The additional line in some cases corresponds to the 
value of the critical visibility in a three setting experiment. 



lead to the violation of local realism with visibility equal 
to 0, as spurious signaling-like effects may lurk in the raw 
data (due to experimental drifts). Fortunately, several 
different approaches for fitting a general physical density 
matrix to the measured data have been developed. There 
are, however, still discussions about the best method. 

We apply the numerical method to the data obtained 
in the four-photon experiment reported in Refs. |40l - |43 |. 
Spontaneous parametric down conversion in combination 
with linear optics was used to observe a variety of four- 
photon entangled states: 



\^{i)) = \i{m 



mf 



1-72 



(10011) + 11100)). 



(9) 



We select three of these states, namely, \4'{^)) = 



|GffZ(7r/4))4, 1^(1/3)) 



1*4), |V'(2/3)) 



\D 



(2)\ 



Moreover, two further states, namely: |G'-ffZ(7r/4))2 and 
\W)z are obtained from the state |'(/'(2/3)) after projective 
measurements. The experimental states were deduced 
from a full quantum state tomography with a maximum- 
likelihood estimation [43| of the most likely density ma- 
trix to have led to the observed data. 

The Steam-roller was used to determine the criti- 
cal parameter v'^^f^ for the tomographically reconstructed 
states (see Tab. IIIII) . If w^J!,^ < 1, the experiment cannot 
be described by any local realistic model. The critical 
parameters v'^ff can be compared with theoretical ones 
''^crit^ obtained for ideal states. The critical visibility in 
the experiment is always larger than the theoretical one. 
This is due to the non-ideal preparation of the states 
in the experiment, which essentially reduces the non- 



classical correlations. Further, the ratios v 



theor 






are 



slightly greater than the experimental fidelity of a given 
state, e.g. for the Dicke state (#57) i'*ri°7^'crft = 0-86, 



No. N State mi 


X 1712 X • • • X 




#49 2 


|GHZ(V4)>2 


2x2 


0.7071 0.7751 


#50 




3x3 


0.7071 0.7737 


#51 3 


\W)3 


2x2x2 


0.6442 0.7534 


#52 




3x3x3 


0.6048 0.7073 


#53 4 


|GHZ(7r/4))4 


2x2x2x2 


0.3536 0.4365 


#54 




3x3x3x3 


0.3408 0.4284 


#55 


l*4> 


2x2x2x2 


0.5303 0.5748 


#56 




3x3x3x3 


0.4823 0.5410 


#57 


\D)f^ 


2x2x2x2 


0.4714 0.5501 



TABLE III. The critical visibilities for some states {theor - 
theoretical and exp - experimental) and number of measure- 
ment settings. 
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Appendix A: Local realistic model for the noisy 

\GHZ{-!v/4))2 state 



whereas F^^f,^ = 0.83. This can be explained by the 
fact, that the noise observed in the experiment contains 
nonclassical correlations [44, EH . 



V. CLOSING REMARKS 

In this paper we used linear programming as a tool 
to study the non-classicality of quantum states. The de- 
tailed comparison of a number of multipartite entangled 
states demonstrates the power of this method even for 
modest computer equipment. Most of the conclusions 
were spelled out in the earlier sections. Here we want 
to stress that the overall message of the obtained data 
is that with more settings per observer one gets stronger 
violations of local realism. The strength of violation also 
increases with the number of qubits (see Fig. [Ij. 

We would like to draw the attention of the reader to 
the peculiarity of three-qubit GHZ states. It is well 
known that the threshold visibility for the GHZ corre- 
lations to violate the Mermin three qubit inequality is 
0.5. This inequality only considers highest order corre- 
lation functions. This threshold is very stubborn. As 
outlined previously, in our numerical approach we used 
not the correlation functions, but the full set of probabil- 
ities, for events of all orders. Still, even if one increases 
the number of settings from two per observer to five the 
stubborn threshold stays put at 0.5, see also [4^. An 
interesting extension of the phenomenon are four-qubit 
GHZ states: if one increases the number of settings from 
two to four for three observers only, with the last one al- 
ways using two settings, the threshold visibility is 0.3536 

as 



for all such cases. The exact value seems to be 



2\/2' 



such is the analytic value for standard Bell inequalities 
(in the two settings per observer scenario). However in 
the 3x3x3x3 case, suddenly the critical visibility drops to 
0.3408. All this points to some so far unexplained resis- 
tance of three-particle correlations, even within a four- 
particle GHZ state, to reveal more non-classicality with 
an increasing number of settings. It should be an exciting 
task to find a reason for that behavior. 



The quantum probability function for the GHZ(7r/4)2 
state has the following form: 

PQAlira,n\0a,(t)a; Ob, (pb) 



= COS (ra 7r/4 -I- 9a) cos (rf,7r/4 -I- 9b) 
+ sin (ra 7r/4 + Oa) sin (rf,7r/4 + 9b) cos ( 



(Al) 



Let us choose some arbitrary set of measurement angles: 
0{ = 0, (/>! = 0; 0? = 7r/2; (/)^ = 7r/4; 0^ = 37r/4 and put 
(|A1|) to ([T]). Then we get the set of constraints, under 
which the visibility function v should be maximized: 



Pa+Pi+ P5 + Pi - 


= q~ 


P2+P3+P7 +P6 


= q+ 


Pi+Pb +P1Z +P12 - 


= q~ 


P& +P7 +P14+P15 - 


= q+ 


P0+P2+P8+ PW - 


= q+ 


Pi +P3+P9+P11 - 


= q~ 


Pi+P6+ Pl2 + Pli - 


= q- 


P5+P7 + Pl3 + Pl5 - 


= q+ 


PS +P9 +P13 +P12 - 


= q+ 


PlQ +Pll +Pli+Pl5 - 


= q~ 


Pa+P2+P6+Pi- 


= q+ 


P1+P3+P7 +P5 - 


= q~ 


PS+PW +Pl2+Pli - 


= q^ 


P9+P11 +P13+P15 - 


= q+ 


Pa+Pi+P9+P8-- 


^q+ 


P2+P3 +Pll +PW - 


= q- 



(A2) 

with q± = (1±w/a/2)/4, Pi < 1 and X)iK = 1- fo'" such 
a problem the numerical procedure gives the highest pos- 
sible visibility Vcrit — 0.707107, for which a local realistic 
model of (|A1[) exists. For such a value oi v = Vcrit, the 
local realistic model has the following form: 



po = 0, pi = 1/8, P2 = 0, P3 = 1/8, 

Pi = 1/8, Pb = 1/8, P6 = 0, P7 = 0, 

Ps = 0, P9 = 0, pio = 1/8, pii = 1/8, 

Pl2 = 1/8, Pl3 = 0, pi4 = 1/8, pi5 = 0. 



(A3) 



8 



Note that in this example, the measurement was cho- 
sen arbitrarily (but in an optimal way). In a standard 
procedure, there is an additional optimization over mea- 
surement angles. 



Appendix B: Local realistic model for a bound 
entangled state 

The local realistic model for the quantum probabili- 
ties for the four-qubit Diir state is presented. The mea- 
surement angles are chosen in the optimal way for the 
GHZ(7r/4)4 state, namely, 9j = 0{i = l,2,3,4;j 
1,2),0} =0,02 ^ ^/2,0i 



7r/4,02 ^ 3^4, 0J 



7r/4, (/)§ = 37r/4, (t)l = 7r/4, ^4 = 37r/4. Any other choice 



of measurements does not have an impact for violation 
of local realism. The model has the following form: 

• PO = Pl5 = P2I = P26 = P38 = Pil = P51 = P60 = 
Pel = P76 = P1I2 = P127 = Pl50 = P153 = Pl65 = 

P170 = P214 = 0.0404029, 

• P85 = P90 = PW2 = PW5 = Pl28 = Pl43 = Pl79 = 
PI88 = P195 = P2OO = P204 = P217 = P229 = P235 = 

P24Q = P248 = P255 = 0.00379126, 

• P197 = P2I9 = P22& = P238 = P245 = 0.0366117, 

• P216 = P233 = 0.0328204. 
All other probabilities vanish. 
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